On the origin of the diffuse extragalactic gamma-ray background radiation 
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We show that inverse Compton scattering of cosmic-microwave-background and starlight photons 
by cosmic-ray electrons in the interstellar and intergalactic space explains well the spectrum and 
intensity of the diffuse gamma-ray background radiation (GBR), which was measured by EGRET 
aboard the Compton Gamma Ray Observatory (CGRO) in directions away from the Galactic disk 
and centre. The Gamma Ray Large Area Space Telescope (GLAST) will be able to separate the 
Galactic foreground from the extragalactic gamma-rays, and to provide stringent tests of the theory. 
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The intensity and spectrum of the diffuse 7 radiation 
was measured by EGRET aboard the Compton Gamma 
Ray Observatory (CGRO). An extragalactic gamma-ray 
'background' radiation (GBR) was inferred PJ from the 
extrapolation of these measurements (in directions away 
from the Galactic disk and center) to zero column den- 
sity, which should eliminate the Galactic contributions of 
bremsstrahlung from cosmic-ray electrons (CREs), and 
7T° production by cosmic ray (CR) nuclei. This GBR 
flux in the observed range of 30 MeV to 120 GeV, shown 
in Fig. [1] is well described by a power-law: 
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The spectral index of the GBR is the same, 2.1 ± 0.03, 
in all sky directions away from the Galactic disk [l|. 
The normalization of the GBR flux in different directions 
was found to be normally distributed around the value 
in Eq. (JTJ. These results were used to argue for a cos- 
mological (extragalactic) origin of the GBR [l[. A large 
number of putative sources have been proposed, from the 
quite conventional to the decisively speculative. Perhaps 
the most conservative hypothesis is that the GBR is the 
sum of 7-ray emissions from unresolved active galactic 
nuclei (AGNs) @. The fact that all AGN detected by 
EGRET are blazars with a power-law 7-ray spectrum 
with an average index 2.15 ± 0.04, compatible with that 
of the GBR, supports this hypothesis [3j , but later studies 
have shown that only < 25% of the GBR can result from 
unresolved AGNs [J]. Geminga-type pulsars, expelled 
into the Galactic halo by asymmetric supernova explo- 
sions, could also be abundant enough to explain the GBR 
||. Other suggestions include cosmic-ray interactions 
in galaxy clusters and groups @, and fossil radiation 
from shock-accelerated CRs during structure formation 
• More exotic hypotheses are a baryon-symmetric Uni- 
verse i, now excluded [§], primordial black-hole evapo- 
ration [lAJ] , supermassive black holes at very high rcdshift 
[llT ] , and the annihilation of dark- matter particles [l2j • 

The EGRET GBR data in directions away from the 
Galactic disk and centre show a significant deviation 




FIG. 1: The GBR spectrum, inferred by EGRET [1]. The 
line is their best power-law fit. 



from isotropy, clearly correlated with the structure of the 
Galaxy and our position relative to its centre [H|. This 
advocates a large Galactic contribution to the GBR. In 
[H[ it was shown that the EGRET GBR could be domi- 
nated by inverse Compton scattering (ICS) of the cosmic 
microwave background radiation (MBR) and starlight by 
Galactic cosmic ray electrons, provided that the Galactic 
cosmic-ray halo is large enough. Indications of a large 
Galactic contribution to the GBR were found by means 
of a wavelet-based "non-parametric" model-independent 
approach Q. Other authors also found that the 
contribution of inverse Compton scattering of starlight 
and microwave background radiation photons by Galac- 
tic CREs is presumably much larger than expected. Ear- 
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lier evidence that ICS by CREs in the Galactic halo con- 
tributes significantly to the GBR at large Galactic lati- 
tudes was reported in [l5| : ~ 30% of the intensity of the 
GBR at large latitudes is correlated to the diffuse Galac- 
tic radio emission at 408 MHz, which is dominated by 
synchrotron radiation from the same CREs that produce 
~ 100 MeV 7-rays by ICS from the Galactic star-light. 

The uniformity of the GBR spectral index over the 
whole sky, despite a large Galactic contribution corre- 
lated with the structure of the galaxy's halo and our po- 
sition relative to its centre, suggests similar origins of the 
Galactic and extragalactic contributions. In this letter 
we argue that ICS of the cosmic microwave background 
radiation by extragalactic CR electrons [IB] , when added 
to the Galactic foreground, explains the EGRET GBR. 
The extragalactic component is calculated directly from 
the CR luminosity of the main putative cosmic accel- 
erators [ifj : supernova explosions and accreting mas- 
sive black holes in active galactic nuclei (AGN). Un- 
like in previous estimates of the contribution from AGN 
0,0,13, which included only the 7-ray emission from 
blazars (AGN with jets and 7 rays beamed in our di- 
rection), we calculate the much larger contribution from 
CREs injected into the intergalactic medium (IGM) in 
arbitrary directions by all AGN jets, and subsequently 
isotropized by the IGM magnetic fields. The Galactic 
component was calculated as in [13[ from our estimated 
Galactic CR luminosity [l6[ and from the locally ob- 
served flux and spectrum of CREs. We show that the 
observed spectrum, intensity and angular dependence of 
the EGRET GBR are correctly predicted. 

The energy spectrum of CREs near Earth [T3], with 
energy E e > 5 GeV, is well described by: 
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The spectral index is predicted by the Cannonball 
(CB) model, wherein CRs are particles of the interstel- 
lar medium accelerated by relativistic "cannonballs" - 
emitted in core-collapse supernova (SN) explosions — to a 
"source" spectrum with a power-law index, (3 S = 13/6 « 
2.17 fl6j . Energy loss by synchrotron emission in mag- 
netic fields and ICS of radiation change f3 s for CREs to 
Pe = Ai+1 ~ 3.17. Radio observations of synchrotron ra- 
diation emitted by CREs in the Galaxy, external galaxies, 
galaxy clusters and AGN, support this predicted univer- 
sal spectrum of high-energy CREs. 

The temperature and mean energy of the MBR are 
T = 2.725 K and e » 2.7 kT w 0.64 meV Starlight 
has e\ ~ 1 eV. Consider the ICS of these radiations by 
CREs. The mean energy of the upscattered photons is: 
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The ICS of the microwave background and starlight pho- 
tons by CREs produces a GBR with a spectrum which is 



a convolution [23j of the CRE spectrum with a thermal 
target spectrum. The result can be approximated by: 
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with E\ obtained from Eq. ([3]) by inverting E n . Intro- 
ducing the electron flux of Eq. into Eq. (|4]) , we obtain: 
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The predicted index agrees with the measured one, 
2.10 ±0.03 [1]. Given Eq. ©, CREs of energy E MBK > 96 
GeV produce the GBR above 30 MeV by ICS of the cur- 
rent (z = 0) MBR; CREs with energy E± > 2.4 GeV suffice 
for ICS on starlight. Let a T «0.65 x 10~ 24 cm -2 be the 
Thomson cross-section and let Ui = rii €4. In our neigh- 
bourhood, U* ~ C/mbr = 0.26 eV cm~ 3 . For electrons 
of energy the radiation cooling times are T ra( j(«) = 
3 to 2 c 3 /(4cr T Ei Ui), so that locally T rad (*) - 6 x 10 s y 
and globally r rad (MBR) - 1.3 x 10 7 /(1 + z) 4 y. These 
numbers are much shorter than a Hubble time. ICS of 
MBR photons dominates the production of the extra- 
galactic GBR, as we argue next. 

Adopt a standard cosmology with H « 70 km s _1 
Mpc^ 1 and (O, Q M , ft A ) = (1, 0.27, 0.73), for which the 
age of the Universe is ~ H^ 1 ~ 14 Gy. For a Galactic 
magnetic field -B^3/iG, Ub ~ Umbk\ synchotron cooling 
and emission are locally relevant [13J ] . In our model, CBs 
transfer their kinetic energy to CRs all along their trajec- 
tories, which extend from the SN-rich inner galaxies to 
their halos and beyond. In galactic halos and galaxy clus- 
ters, B<3fiG, and in the IGM, B ~ 50 nG Q. In both 
places starlight is irrelevant, and ICS of the MBR whose 
energy density increases with z like (1 + z) 4 dominates 
over synchrotron losses on the magnetic field. Thus, we 
calculate the intensity of the extragalactic GBR from the 
conclusion that the kinetic energy of CREs in the Uni- 
verse with a lifetime shorter than the Hubble time has 
been converted by ICS of the MBR to 7-rays with the 
predicted spectrum of Eq. |5} . 

In the CB model, the main accelerators of high- 
energy CREs are the relativistic jets of supernovae (SNe) 
and AGNs [l6j]. The SN rate is proportional to the 
star-formation rate i? SN (z) cx R SF (z), with i? SN (0) w 
10~ 4 Mpc -3 yr _1 The observations are well rep- 

resented by i? SF (z)/R SF (0) ~ (1 + z) 4 for z < 1 and 
R SF (z) = i? SF (l) for z > 1.2 |3| • Let £ t «2x 10 51 
erg be the mean energy release in CRs per SN [l6[ and 
let / e be the fraction of the luminosity in CREs out of 
the total L CR in CRs. The CB model does not predict 
f e , we assume that it is equal to the ratio of the Milky 
Way's luminosity in CREs to its total luminosity in CRs: 

_ C <±F e [MW] EdE 1 f dF CR [MW] EdE _ 1 
/e ~J dE ^~7/ J dE ~ 40 ( } 

where r csc « 2 x 10 s (E /GeV) -0 ' 6 yr is the mean escape 
time of CR protons and electrons from the Galaxy and 
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its halo by diffusion in its magnetic field [16[ . In the CB 
model the volumes occupied by electron and proton CRs 
are similar, because CBs generate CRs all along their 
trajectories, which constitute a dense mesh in the Galaxy 
and its halo. The integrals in Eq. ([6]) extend from a lower 
fixed Lorentz factor, which drops from the ratio f e , since 
the integrands are source spectra, identical in the CB 
model for electrons and protons. 

The energy of CRE made by SN jets is converted, 
above E e ~ 100 MeV, to photon energy. Their contri- 
bution to the GBR spectrum satisfies: 
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where L e = f e L CR ~/ e 1/6 ~ / e i? SN ^ fc 7, 
luminosity density of CREs with Lorentz factor above 
7 C = y/3 E c /4 eo- Hence, Eq. ||7J) yields for the contribu- 
tion to the GBR from extragalactic SNe: 
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Powered by mass accretion onto massive black holes, 
AGNs eject powerful relativistic jets whose kinetic en- 
ergy is transferred mainly to CRs. The kinetic power 
of these jets has been estimated from their radio lobes, 
assuming equipartition between CR- and magnetic field 
energies and an energy ratio f e similar to that observed 
in our Galaxy. It was estimated [2l[ that AGNs with a 
central black hole of M ~ 10 8 M inject « 10 61 " 62 erg 
into the intergalactic space, mostly during their ~ 10 s y 
bright phase around redshift z = 2.5. In search for an up- 
per bound, we assume that the kinetic energy release in 
relativistic jets is the maximal energy release from mass 
accretion onto a Kerr black hole (~ 42% of its mass), 
and that this energy is cquipartitioncd between magnetic 
fields and cosmic rays [16[ with a fraction f e of the CR en- 
ergy carried by electrons. These CREs also cool rapidly 
by ICS of the MBR. The energy of CREs whose radiative 
cooling rate T ra d(z) is larger than the cosmic expansion 
rate, H(z), is converted to 7-rays. Their energy is red- 
shifted by 1 + z by the cosmic expansion. Using a black 
hole density, p^Az = 0)~2x 10 5 M Q Mpc -3 in the cur- 
rent Universe |22j |. and the CB-model injection spectral 
index [l6j , our estimated contribution from AGNs to the 
extragalactic GBR flux, 
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This upper limit is smaller than the SN result of Eq. (jHJl . 

The GBR contains a considerable Galactic foreground 
due to ICS of MBR, starlight and sunlight photons by 
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FIG. 2: The flux of GBR photons above 100 MeV: comparison 
between EGRET data and our model for h e = 8 kpc, p c = 35 
kpc, as functions of longitude I at fixed latitudes b. The 
shaded domain is EGRET's mask. Notice that the vertical 
scales do not start at zero. 

Galactic CREs 16]. The convolution of a CRE power- 
law spectrum with a photon thermal distribution [231 ] can 
be approximated very simply [l3j |. Using the index i to 
label the MBR, starlight and sunlight fluxes, we have: 
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where Ni(b, I) is the column density of the radiation field 
weighted by the distribution of CREs in the direction 
(6, I), and E\ is given in Eq. ((4]). The distribution of 
the non-solar starlight is approximated as oc 1/r 2 , with 
r the distance to the Galactic centre, and the CREs are 
assumed to be distributed as a Gaussian "CR halo" [HI . 
Naturally, the results depend crucially on the size and 
shape of this halo. In this note we use our updated esti- 
mate of the CR halo [l6j]: a Gaussian distribution with 
a scale l eng th of p e = 35 kpc in the Galactic disk, as we 
used in [131 ] , but a scale height of h e = 8 kpc perpendic- 
ular to the disk [l6[ instead of the h e — 20 kpc used in 
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[T3 |. The justification for this change is as follows: The 
radio emission of "edge-on" galaxies -interpreted as syn- 
chrotron radiation by electrons on their magnetic fields- 
offers direct observational evidence for CREs well above 
galactic disks (e.g. [25]). For the particularly well ob- 
served case of NGC 5755, the exponential scale height 
of the synchrotron radiation is 0(4) kpc. If the CRs 
and the magnetic field energy are in equipartition, they 
should have similar distributions, and the Gaussian scale 
height h e of the electrons ought to be roughly twice that 
of the synchrotron intensity, which reflects the convo- 
lution of the electron- and magnetic-field distributions. 
The inferred value is h e ~ 8 kpc. The corresponding 
volume of the Galactic CR halo is V CR = (tt) 3 ^ 2 pi h e — 
1.6 x 10 69 cm 3 . The SN rate in the Galaxy is i? SN [MW] ~ 
2 per century, and its predicted total luminosity in CRs 
is L CR w Ek i? SN [MW] « 4 x lO^ergy" 1 . The CR con- 
finement volume must obey the constraint: 

4tt f dE dF v , s 

i «~"«x T j—B-g. (12, 

Our estimated r osc and the observed (or fitted) spectrum 
of CRs [H [SI yield the expected V CR w 1.6 x 10 69 cm 3 . 
The volume inferred from a leaky-box model fit to the 



Galactic GBR (2?J is smaller by a factor w 2.5 than our 
estimate, reflecting the shorter confinement time of CRs 
estimated in leaky-box models from the abundance of 
unstable CRs (28|, and the higher contribution assumed 
in [27J for the extragalactic GBR. 

In Fig. 2 we compare the observed GBR with our 
predictions, as functions of Galactic coordinates. The 
prediction is a sum of a (b, Z)-dependent Galactic fore- 
ground produced by ICS of the MBR, starlight and sun- 
light, and a uniform extragalactic GBR. The result has 
X 2 /dof = 0.85, a vast improvement over the constant 
GBR fit by EGRET, for which x 2 /dof = 2.6. The ratios of 
/-integrated extragalactic to galactic fluxes are ~0.5, 0.9, 
1.5, for \b\ = 20°, 45°, 75°. The 'foreground' component 
of the 7 'background' is ~50% of the total radiation. 

We conclude that the GBR can be explained by stan- 
dard physics, namely, ICS of MBR and starlight by CREs 
from the two main CR sources in the universe: SNe and 
AGNs. At E 1 > 100 GeV, most of the extragalactic GBR 
is absorbed by pair production on the CMB [2^] and the 
diffuse GBR reduces to the Galactic foreground. This 
suppression, and a decisive determination of the angular 
dependence in Fig. 2, should be observable by GLAST. 
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